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We present a study of D — > /T, Iv semileptonic decays on the lattice which employs the HISQ ac- 
tion for both the charm and the light quarks. We work with MILC unquenched A'/ = 2 + 1 lattices 
and determine the scalar form factor fo{q^)- This form factor is obtained from a scalar current 
matrix element that does not require any operator matching. We find /o'^^(O) = /+^^(0) = 
0.747(19) in the chiral plus continuum limit and hereby improve the theory error on this quan- 
tity by a factor of ^4 compared to previous lattice determinations. Combining the new the- 
ory result with recent experimental measurements of the product /+^^(0) * jVci l from BaBar 
and CLEO-c leads to a very precise direct determination of the CKM matrix element |Ve.v|, 
I Veil = 0.961 (11) (24), where the first error comes from experiment and the second is the lattice 
QCD theory error 
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1. Introduction 

From a study of D ^ K, Iv semileptonic decays, one can calculate the form factor /+ {q^ = 0). 
One can also determine the CKM matrix element, |Vc^|, by combining theory and experimental 
inputs. We continue to work on the D semileptonic decay project that was presented at the Lattice 
2009 conference In this article, we present a brief summary of our recent results for the Dto K 
semileptonic decays, which are already published in Ref. [^]. So, for more detail, please see the 
publication. 

For this project, we use A^y = 2 + 1 asqtad MILC gauge configurations with two lattice spac- 
ings, a ~ 0. 12fm "coarse" and a ~ 0.09fm "fine" ensembles. We apply the HISQ action for both the 
charm and light valence quarks. For better statistics, we employ random wall sources. We develop 
a new extrapolation method to go to the continuum and chiral limit, the so called "simultaneous 
modified z-expansion extrapolation," which allows us to extrapolate the form factors for the en- 
tire q'^ range. This method does not have the expansion problem which normal chiral perturbation 
theory would have at large Ek- 

To study the process D ^ K,lv one needs to evaluate the matrix element of the charged elec- 
troweak current between the D and the K meson states, {K\{V^ —A^)\D). Only the vector current 
contributes to the pseudoscalar-to-pseudoscalar amplitude and the matrix element can be writ- 
ten in terms of two form factors f+{q^) and f{){q^), where q^ = p^ — is the four-momentum of 
the emitted W-boson. 



{K\V^\D)=f^-^^{q' ^ 

mI-m} 



(1.1) 



with = Ysj^^c. As described below, we find it useful to consider also the matrix element of 
the scalar current S = Ys^o 

{K\S\D) = ^^i^/r'iq')- (1-2) 

moc - mos 

In continuum QCD one has the PCVC (partially conserved vector current) relation and the vector 
and scalar currents obey, 

q^{V;,""'-) = {m,,-mo.){S'''"'-). (1.3) 
In fact PCVC is the reason why the same form factor fg^^iq^) appears in eqs. (|n|) and ([0|). On 



the lattice it is often much more convenient to simulate with vector currents YqiY^^q2 that are not 
exactly conserved at finite lattice spacings even for Ql = Q2. Such non-exactly-conserved currents 
need to be renormalized and acquire Z-factors. We are able to carry out fully nonperturbative 
renormalization of the lattice vector current by imposing PCVC. In the D meson rest frame the 
condition becomes, 

{Md - Ek) (C- )Zr + pK- {V'""- )Z, = {moc - mo.) {S'^"' ) . (1.4) 

We have checked the feasibility of this renormalization scheme and extracted preliminary Zt 
and Zs values for the test case of Dg — )• rj^Jv in Ref.[jl|]. However, here we focus on the form factor 
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f+{q^) just at = 0, since this is all that is needed to extract \Vcs\- We do this by exploiting the 
kinematic identity /+(0) = /o(0), and concentrating on determining the scalar form factor fo{q^) 
as accurately as possible. The best way to proceed is to evaluate the hadronic matrix element of the 
scalar current rather than of the vector current. From eq.( |l .2| ) one then has, 

' Ml- Ml 

The numerator on the right-hand-side is a renormalization group invariant combination. This is 
true even in our lattice formulation, because we use the same relativistic action for both the heavy 



and the light valence quarks. Moreover, eq. (1.5) allows a lattice determination of fo{q ) and hence 



also of /+(0) = /o(0) without any need for operator matching. Using eq. (1.5) and going to the 



continuum limit is straightforward, because our action is so highly improved even for heavy quarks. 

2. Simultaneous modified z-expansion extrapolation 

The continuum z-expansion method is a well known model-independent parameterization 
method for semileptonic decay form factors. One can write the form factor as, 

where P{q^) and ^o{q^,to) are given functions from analyticity properties of the form factors. 

The z-expansion method works well for individual ensembles, however we like to modifying 
the fit ansatz to enable extrapolation to the physical limit. All kinematic properties that depend 
on q^ are absorbed by P,'i>o, and z. A natural way to distinguish between ensembles is to let 
Qk ^ aji* Du, where D^- contains the light quark mass and lattice spacing dependence as shown 
below with kmax = 2. 

/o(^^) = («o£>o + aiDiz + a2£>2Z^) (2.2) 



x{\ +bi{aEKf + b2{aEK^^^ 



where. 



Di = \ + c\xi + c^dxs + c\xilog{xi) + di{amc) (2.3) 
+ei{am,f + fi(^5Ml + 5Ml 



In eq. |2.3| , we put typical analytic terms for light valence (x/ and 5xs terms) and sea quark mass 
{5M;i and 5Mk terms) dependence. For the chiral logs, we only include up/down quark contri- 
butions. The strange quark chiral logs are close to a constant that can be absorbed into the a,'s. 
There are two distinct sources of lattice spacing dependence. {arricY and {anicY terms are due to 
the heavy quark discretization error, and {qEkY and {aEKY terms are introduced to estimate the 
discretization errors due to finite momentum. Since we want the a,D/ to be independent of the 
momentum, the aEx terms are placed separately outside the z-expansion. We include lattice spac- 
ing dependent terms up to fourth power, however we tested with even higher terms and confirmed 
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Figure 1: Chiral/continuum extrapolation of fo{q^) versus E-^ from the modified z-expansion ansatz. The 
data points are coarse (left) and fine (right) lattice points. Three individual curves and the extrapolated band 
are from a fit to all five ensembles. 



that the higher terms are negligible. We have carried out simultaneous fits to all the data using the 
above ansatz and find that very good fits are possible. Fig. |l] shows the resulting fit curves for each 
ensemble and the chiral/continuum extrapolated curve with its error band for fo{q^) versus (we 
show separately the coarse and fine ensembles in order to avoid too much clutter). On the left panel 
of Fig. ^ we show /o {q^ = 0) for the five ensembles and in the physical limit. One sees that within 
errors this quantity shows little light quark mass dependence and a ~ 1.3% lattice spacing depen- 
dence. We also test the chiral/continuum extrapolation with partially quenched chiral perturbation 
theory (PQChPT). This traditional method gives results in very good agreement with the modified 
z-expansion extrapolation method (see the right panel of Fig. 

3. /+(0), \ Vcs\, and unitarity tests 

3.1 /+(0) = /o(0) 

From the simultaneous modified z-expansion extrapolation method, we find /+(0) = 0.748 ± 
0.019 in the physical limit for D" R-lv, and /+(0) = 0.746 ±0.019 for D+ Z"/v. We take 
an average over these two channels and our final result in the physical limit becomes, 

/{^^^(O) =0.747 ±0.011 ±0.015. (3.1) 

The first error comes from statistics and the second error represents systematic errors. Table |] 
summarizes the error budget. One sees that the largest contributions to the total error come from 
statistics followed by {anic) and {aEx) extrapolation errors. 

In order to calculate the form factor, we have to put in meson masses from experiment and also 
from our lattice simulations. For example, we need experimental D, K, and TT meson masses to get 
the form factor at the physical limit, and Ek, D, and K meson masses from the lattice calculations 
are used to fit at non-zero lattice spacing. In Table |l|, "Input meson mass" refers to errors induced 
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Figure 2: (left) /o at g'^ = for the five ensembles and in the physical limit, (right) Comparisons of fQ{q^) 
in the physical limit from the z-expansion and the ChPT extrapolations. 



Type 


Error 


Statistical 


1.5 % 


Lattice scale (ri and ri/a) 


0.2% 


Input meson mass 


0.1 % 


Light quark dependence 


0.6% 


Strange quark dependence 


0.7 % 


Sea quark dependence 


0.4% 


ante extrapolation 


1.4% 


uEk extrapolation 


1.0% 


Finite volume 


0.01 % 


Charm quark tuning 


0.05 % 


Total 


2.5% 



Table 1: Total error budget. 



from these input meson masses. In the fit ansatz, eq. 23, there are light quark (c\ and C3), strange 
quark {c'2), and sea quark dependent terms {ft). Each systematic error due to these terms is shown 
on the fourth to sixth line in the table. Lattice spacing dependence errors are estimated separately 
for {arricY and {uEk)-' type contributions. 

In the fit ansatz, xilog{xi) is the most infrared sensitive term. We calculate the pion-tadpole 
loop integral both at finite volume and at infinite volume and compare these to estimate the finite 
volume effects. For the charm quark mass tuning error, we calculate the form factor with a different 
charm quark mass, anic = 0.629, on the C3 ensemble, and compare with the result with the tuned 
ante = 0.6235. 

In their papers both BaBar [Q] and CLEO-c have converted their measurements of /+(0) * 
\Vcs\ into results for /+(0) using values for \Vcs\ fixed by CKM unitarity. For this CLEO-c uses 
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HPQCD (2010) 

Other Iheury 
Fermilab/MILC (2005) h 

Sura Rules (2009) H 



Experimeiil + CKM Unitarity 
CLEO-c (2009) 
BaBar (2007+update) 



B 2.5 % error - 
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f,(q'=0) 



-| ^ 1 ^ 1 ^ 1 ^ 1 r 

HPQCD (2010) 
IV 1 = 0.961 (11)(24) 



-B- 



PDG 2010: direct estimation 
Semi-leptonic decay: \VJ = 0.98 (1)(10) • 

Leptonic decay: IV^J = 1 .030 (38) i-qh 



Average: IV^,^I = 1.023 (36) 

PDG 2010: CKM Unitarity 
IV 1 = 0.97345 (16) 



J , I , L 



i-en 
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IV I 



Figure 3: (left) Comparisons of /()(^^ = 0) with other calculations and experiments, (right) Comparisons of 
our new \Vcs\ with values in the PDG [H]- 



the 2008 PDG CKM unitarity value of |Vc.v| = 0.97334(23) and obtains /f^'^(O) = 0.739(9) and 
BaBar uses = 0.9729(3) leading to /+(0) = 0.737(10). On the left panel of Fig. | we plot 
our result, eq.(3.1), together with earlier theory results from the lattice ^ and from a recent sum 



rules calculation and with the BaBar and CLEO-c numbers. One sees the very welcome reduction 
in theory errors which are now small enough so that the agreement between theory and experiment 
already provides a nontrivial indirect test of CKM unitarity. 

3.2 Direct Determination of \Vcs \ and unitarity tests 

As experimental input we take /+(0) * |V„| = 0.719(8) from CLEO-c [|] and /+(0) * |Vc.v| = 
0.717(10) from BaBar [Q]. For the latter we have multiplied BaBar's quoted /+(0) with their 
quoted CKM unitarity value for |V„j. Averaging between the two experiments we use /+(0) * 



|V„| = 0.718(8) together with eq.(^ to extract IV^I- One finds, 



|y„ I = 0.961 ± 0.01 1 ± 0.024, (3.2) 



in good agreement with the CKM unitarity value of 0.97345(16) The first error in ( [3.2| ) is from 
experiment and the second from the lattice calculation of this article. This is a very precise direct 
determination of \Vcs\, made possible by the many advances in lattice QCD that are described in 
this article together with the tremendous progress in recent experimental studies of D semileptonic 
decays [Q, |5|]. On the right panel of Fig. ^ we plot several previous direct determinations of \Vcs\ 
from the 2010 PDG ||] together with (O) and the CKM unitarity value. 



Using the new value of [V^l, eq.(3^), and the current PDG values \Vcd\ = 0.230(11) and 
iVcfal = 0.0406(13) one finds. 



\Vcd\^ + \Vcs\^ + \Vch\^ = 0.978(50) (3.3) 
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for the 2nd row. And similarly for the 2nd column, with \Vus\ = 0.2252(9) and \Vts = 0.0387(21) 
one gets, 

|K.vP + |V„.|2 + |y,,|2 = 0.976(50). (3.4) 

4. Discussion 

We have carried out a successful calculation for D ^ K,lv semileptonic decay form factors 
using the HISQ action for both the charm and light quarks with Nf = 2 + 1 asqtad MILC gauge 
configurations. The total error for /+(0) is estimated here to be 2.5%. This is a factor of four 
times smaller than in the previous lattice calculation of Ref. [Q]. This was achievable because of 
applying several new methods and techniques. We employ the HISQ action for both charm and 
light quark actions and a scalar current rather than the traditional vector current. Because of these 
new methods, we obtain results with smaller discretization errors and no operator matching. We 
also developed the modified z-expansion extrapolation method, which is crucial to decrease errors 
due to the discretization, chiral / continuum extrapolation and parameterization of the form factor. 
In order to decrease statistical errors, we apply random-wall sources and perform simultaneous fits 
with multiple correlators and T's. If we compare with the error budget of Ref. 0, then we see the 
statistical errors reduced from 3% to 1.5% and the extrapolation and parameterization errors from 
3% to 1.5% as well. The biggest improvement is in the discretization errors. The total discretization 
errors have now been reduced from 9% to 2%. We note that the concept of the discretization 
errors is different in Ref. compared to ours. In Ref. [^, they estimate the discretization errors 
by power counting, since they calculate at only one lattice spacing. Here, however, we actually 
perform continuum extrapolations with correction terms for the discretization effects. As a result, 
we do not have discretization errors per se, but instead extrapolation errors due to higher order 
correction terms. 

Again, this is a short version of Ref. [^. For more detail and full discussion, please see the 
publication. 
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